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Abstract 

Recently, a class of gravitational backgrounds in 3 + 1 dimensions have 
been proposed as holographic duals to a Lifshitz theory describing critical 
phenomena in 2 + 1 dimensions with critical exponent z >!. We numerically 
explore black holes in these backgrounds for a range of values of z. We find 
drastically different behavior for z > 2 and z < 2. We find that for z > 2 
{z < 2) the Lifshitz fixed point is repulsive (attractive) when going to larger 
radial parameter r. For the repulsive z > 2 backgrounds, we find a continuous 
family of black holes satisfying a finite energy condition. However, for z < 2 
we find that the finite energy condition is more restrictive, and we expect only 
a discrete set of black hole solutions, unless some unexpected cancellations 
occur. For all black holes, we plot temperature T as a function of horizon 
radius tq. For z ^ 1.761 we find that this curve develops a negative slope for 
certain values of rg possibly indicating a thermodynamic instability. 



1 Introduction 



Ever since the Maldacena conjecture [1], holography has become an important tech- 
nique for studying strongly coupled systems (for a review, see [2]). Traditionally, 
much of the work has been devoted to the study of relativistic field theories in 3 + 1 
dimensions with superconformal symmetry [3]. However, holographic methods have 
also proven useful in considerably less symmetric situations [4]. It has also proven 
useful to construct toy models [5] imitating the string theory backgrounds, with the 
surprising result that much of the physics is capturedQ 

From these 3 + 1 dimensional setups, it is known that when describing finite 
temperature field theories one must consider the black holes in the dual geometry. 
The canonical example is that the thermodynamic properties of AdS-Schwarzschild 
black holes [6] match thermal properties of A/" = 4 SYM theories [7] , exhibiting a host 
of interesting effects. The black hole backgrounds themselves have led to interesting 
conjectures about the nature of strongly coupled plasmas [8]. There have also been 
recent discoveries relating higher derivative corrections in gravitational actions to 
unitarity in the field theory [9]. 

Given this, one may be curious about what other types of field theoretic systems 
can be modeled in terms of a gravitational theory. Other than the 3 + 1 dimensional 
setups mentioned above, one can also use holography to study 2 + 1 dimensional field 
theories relevant for condensed matter systems. Recently, much effort has gone into 
describing quantum critical behavior for these theories using holographic techniques, 
for a review, see [10]. Although not the only type of system one could study, quantum 
critical systems exhibit a scaling symmetry 



similar to the scaling invariance of pure AdS {z = 1) in the Poincare patch. From a 
holographic standpoint, this suggests the form of the spacetime metric 



where the above scaling is realized as an isometry of the metric along with r r. 
Other metrics exist with the above scaling symmetry, but also with an added Galilean 
boost symmetry [11-14]. 

Often, a good place to begin studying any system is to write down a toy or 
"phenomenological" model [11,12] to study generic properties (for thermal versions 

^although one prefers string theory "brane setups." where one knows how to describe the weakly 
coupled degrees of freedom. 
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of these models, see [15]). One may then consider possible embeddings into a more 
fundamental theory [13, 15, 16], such as a string theory, where more information 
is known about the weakly coupled physics. One could also consider gravitational 
theories which inherently have some nonrelativistic scaling built in [17]. However, 
here we will be content to study the relatively simple model considered in [18] where 
Kachru, Liu and Mulligan consider an action which admits a solution with metric ([2]) 
(earlier studies of these metrics in a "brane world" scenario appear in [19], and further 
investigated in [20]). We leave the possible embedding of this model into string theory 
for future work, although other kinds of generalizations appear in [21-23]. It should 
also be noted that a related system [16] with anisotropic space scaling has been 
constructed (in 4+1 dimensions) complete with analytic black brane solutions: this 
may serve as a template for embedding these types of theories into string theory. 

It is the current aim of this work to numerically construct the black hole back- 
grounds that asymptote to the metric ([2]) given the equations of motion for the 
model [18]. Our work can be thought of as complementary to the work of [24], [25] 
where they consider the case z = 2. 

We now turn to the model of [18]. The action that we consider is 



S = I d^x^g (R-2K- \t^,T^^ - jA.A'' ) (3) 



where JF = dA. Up to a Legendre transform, this is equivalent to the action given 
in [18]. This can be seen directly from the equations of motion given in [18] with the 
identification that = ^A where F3 is the Kachru et. al. three-form, and A is 
our one-form given above. Further, we will parameterize the constants in the action 
above as 

V2Z , IZ^ + Z + A 
c=^—, A = . (4) 

Note that this allows for arbitrary negative A, and allows for the ratio of physical 
constants to be in the regime | < ^ < cxd. The specific choice ([4]) will become 
convenient shortly. 

There are several known solutions to the above action. First, we consider the 
"black brane" and "black hole" solutions in pure AdS, 

ds^ = (Zpi (-(a + r^f{r))dt^ + r\dxl + {l-acos\x^))dxl)+ '^''^ 



^ = 0. (5) 
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Figure 1: The various fixed points for a given action parameterized by Z and L. 
Tlie (dasfi-dot) green curve is z = ^, the (long-daslied) red sloped line is z = Z and 
tfie (sliort-daslied) blue flat line is z = 1. We choose to plot only Z > 2 because 
these specify different actions that have different ratios We will find that for 
the critical exponent z < 2 (the dash-dot green curve) that the Lifshitz fixed point 
is "attractive" when going to larger r and that for z > 2 (the red sloped line) it is 
repulsive. Further, for Z > 4 the solution z = 4/Z < 1 is imaginary, and so we 
exclude it. 



Here, and throughout, we will be considering two cases simultaneously: a = is 
the "black brane" case and cr = 1 is the "black hole" solution. We will color code 
a and all terms multiplied by this factor in red, so it is easy to see how the various 
equations are modified (however, a will always appear, so that the color coding is 
redundant). 

There is, of course, also the solution discussed in [18], given by 

els'" = L^(^-r^^dt^ + r'idxi + dxl) + ^'^ 
l2z(z - 1) , 

A = L'-^ \, U t (6) 

with the identification z = Z and L = L. One may of course invert these equations 
and find 2;(A, c), L(A, c), perhaps more intuitively regarding the parameters of the 
solution as depending on the parameters of the theory. However, here we will only 
be concerned with the restricted set of theories given by (j4]). The above background 
will serve as the asymptotic form of all our black hole/brane spacetimes as r — >^ cx). 
We also note that there may be multiple Lifshitz fixed points given the same 
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action. The most efficient way to find these is to write down two solutions of the 
form ([6]) and requiring that even though the values z, L and z', V are different, the 
parameters of the action c^{Z, L)\^^^ = c^{Z, L)\^^^, i^^, and A{Z, L)\^^^ = 

A{Z, I/) 1 2^^, i^j^, remain unchanged. A second solution is then found, given by 

Hence, the solution 

ds' = L'^(^-r'''dt^ + r\dxl + dxl) + ^'^ 

A ^ (8) 

is also a solution with the parameters given in (jlj). We can see that z,z' > 1 for 
the solutions ([6]) and ([8]) to be real. Therefore, there are two solutions for the 
range 1 < Z < 4 because both Z and 4/Z are bigger than 1. This gives us the 
following picture: we parameterize the constants of the action as (jlj), and so Z and 
L are our free parameters. There exist solutions of the form with the assignment 
z = Z,L = Lorz = ^ and L = ^L. The fixed point of this transformation between 
solutions is Z = 2, and so to avoid double counting solutions, we consider Z > 2, 
which is also the range denoting different actions with different ratios Also note 
that pure AdS is always a solution to the above equations as well (in the global patch 
when using the sphere) and so 2; = 1 with Z unconstrained is also a good fixed point. 
Hence, we obtain a plot of fixed points for various Z in figure [H 

The remainder of the paper is dedicated to the analysis of the equations of motion 
resulting from ([3]). We organize this as follows: in the next subsection we write down 
an Ansatz and the resulting differential equations coming from the action ([3]). In 
section 2, we explore series expansions near a regular horizon (r = tq), and near the 
asymptotic region, r = 00. Further, we discuss the finite energy condition, defined 
as the on shell Hamiltonian for our action [26]. In section 3 we combine the results 
of section 2, and numerically integrate the equations of motion. We end section 3 
with a discussion of the results, some open questions, and possible future directions. 
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1.1 The Ansatz and reduced action. 



As explained above, we will be concerned with constructing the black brane/hole 
type solutions that asymptote to Qj. We do so by considering the Ansatz 



Again, in the above, we have introduced two distinct cases: when a = 1 the two 
dimensional metric is that of the unit 5*^ and when a = it is simply flat two 
dimensional space. To obtain the equations of motion, we will first reduce the action 
to one dimension, plugging in the above Ansatz into action Q. After integration by 
parts, this gives a reduced action 



where we use the abbreviation that d = This action reproduces all Einstein- 
Maxwell equations, if one includes the equation of motion for C(r). C{r) acts as a 
Lagrange multiplier imposing the "zero Hamiltonian" condition. The above action 
contains 3 dynamical fields, and one constraint equation. We expect, therefore, 
6 — 1 = 5 integration constants associated with solutions to the equations near 
a generic point. However, we expect to be able to remove one of the remaining 
constants via coordinate transformations, and so we expect 6 — 2 = 4 constants. 
Using coordinate transformations to fix some relation we will refer to as gauge fixing. 
The above action gives the following equations of motion 



-exp{2A{r))df + exp{2B{r)){{dx'y + (1 - a cos\x')){dxY) 
+ exp{2C{r))dr^ 

A = exp G{r)dt. 



(9) 



(10) 




(11) 



(12) 
(13) 



0. (14) 
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where the last equation is the Hamiltonian constraint. We will be concerned with 
turning on a "blackening" factor for those backgrounds found in [18]; these back- 
grounds all have c 7^ 0. However, in the above we can see why the c = (massless) 
case is special. In such a case, one can shift A ^ A + 6 and C ^ C — 6 with 6 being 
an arbitrary function of r. Doing so leaves the equation of motion for G given by 
(fT3l) unaffected (the Laplacian part remains unchanged). In such a situation, G may 
be left alone, and still satisfy its equation of motion. The mass parameter c changes 
this because of the relative minus sign in front of C(r) in the exponential. Hence, 
we expect the massive vector field above to change in a "blackened" background. 

Before trying to find solutions to the above equations, we present a first integral 
for the (7 = version of the above equations of motion. One may infer the existence 
of this first integral because a = leaves only 2 potential terms in the effective 
Lagrangian. However, because there are 3 fields some linear combination of them 
does not couple to the potential. One may also see this as the Noether charge 
associated with the shift 



/Mr)] 
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with 5 a constant. The above represents a diffeomorphism which preserves the volume 
element dtdx^dx^. This is why it is inherited as a Noether symmetry in the reduced 
Lagrangian. 

Rather than writing the conserved quantity, we write out the differential equation 
that may be directly integrated in the cr = case: 

d{2e^A+2B~C)Qj^ _ 2g{A+2B-C)^^ _ g(-A+2B.C+2G)^^) ^ a2e^^^''^ = 0. (16) 

As an example, (for cr = 0) we integrate once, and plug in the black brane in AdS 
case to find 

SO the above equation is indeed satisfied. 

Given the above considerations, we may pick as a complete set of differential 
equations ([HD, (fT6l) and one of the three second order differential equations (fTTj) or 
(fT2l) or (fT3l) . In fact, one can gauge fix one field to a known function (e.g. = r). 
This will transform one of the three second order differential equations into a first 
order differential equation. In the a = case, one is therefore left solving 3 first order 
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differential equations (after integrating ( fT6l) once). While this seems suggestive, we 
have been unable to use this fact to solve the equations exactly. 

In general, we will want to explore solutions to the equations of motion that 
asymptote correctly to ([6]). For this reason, we take the following definitions 

A{r) = In(rU) + Ai{r), B{r) = In(rL) + Bi{r) 

C{r) = In (^^^ + Ci(r), G(r) = In (^^^j^^i^S^ + Gi(r) (18) 

As mentioned above, the 3 second order differential equations may be reduced to 
two first order and one second order differential equations by a gauge choice, and 
the presence of a conserved Hamiltonian. In the black brane {a = 0) case, one can 
further reduce to 3 first order equations due to the scaling symmetry, however we will 
deal with the a = 1 case here. While we could make many different gauge choices, 
what will be important for us later is the gauge choice 

5i(r) = 0. (19) 

In such a gauge choice, the second order differential equation for B becomes a first 
order equation. In addition, we take the Hamiltonian, and the second order equation 
for e*^. Using the first two of these three equations, we may eliminate de"^"^ and de'"^ 
in the third, which we do. Further, we isolate de^^ and (9e*^i in their own equations. 
Doing so, we find the three equations 



de^^ + - 

4 ze^i ■ 2 e^i 

1 2^62^1 (1 + 2z) - ze^^'e^^' {z^ + z + 4) + z {z - 2e^^') e^^^ {z - 1] 
4 zre"^^ 

1 e^'e^^' 



^2-^ = (20) 



_ i re^i (ge^^)'(z-l) _ 1 e^^e^^ge^Hz - 1) 

^ ' 4 2e2^i 2 e2^i 

1 e<^i (eze^^i - ze^^^e^'^i {z^ + z + 4) + z [z + 2e^^^) e^^^ {z - 1)) 

4 zre^^'^ 

le^^i , , 

+^27^ = (21) 



{z - 1)96^1 (2(2 + l)e2^i - e2^ie2Gi(2 - 1)) 



z{z - l)e^i (-e2^ie2Gi(z - 1) + e^^'{z + 1) - 2e^^^e^^^) 



(22) 
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where we have exphcitly left the z — 1 multiplying the second order equation, to 
stress that this equation does not need to be solved in the z = 1 case. Further, one 
may combine the first two lines into the following equation 

rd (e^^e^i) + e^^e^^ (^1 - (j-^^ j - 1) = 0. (23) 

This merely serves as a check: in the case z = 1 the black brane/black hole solutions 
have e^^e'^^ = 1, which solves the above equation. We will use fl23|) to find an exact 
solution for z = 4 in appendix |X1 

The initial conditions needed to solve the above equations numerically are initial 
values for e^^ , e*"^ , e*^^ and an initial value for de*-^^ . 



2 Analytic explorations. 

2.1 The perturbed solution near the horizon 

We begin by first exploring the solution near the horizon. We require that e'^^ goes to 
zero linearly, e^*" has a simple pole, and goes to zero linearly to make the flux dA 
go to a constant (in a local frame or not). Further, we take the gauge B{r) = In(Lr) 
for this section. We expand 

- ro)^ + aoai{r - Tq)^ H j j , B{r) = In(rL) 

-ro)"5+ci(r-ro)5 + --- n , (24) 




^^^^^2 {(^ogo{r - ro) + aogi{r - r^f H ) 

Note that by scaling time we can adjust the constant by an overall multiplicative 
factor (note the use of Cq in the expansion of G{r) as well, as multiplies dt for 
the one- form A). We will need to use this to fix the asymptotic value of A{r) to be 
exactly In(r^L) with no multiplicative factor inside the log. 

We plug this expansion into the equations of motion arising from (fTOj) . and solve 
for the various coefficients. We find a constraint on the 0th order constants: as 
expected not all boundary conditions are allowed. We solve for Cq in terms of the 



8 



other go and ro, and findS 



^ ^ v/(2z + ^gro(^-"I))r| 



ai 



The higher order coefficients are 

rg(^ - 1)^(7^ + 2ro{z - l)(rg(-z' + 2^ + 4) + a)g^ - Azjrlzjz^ + ^ + 4) + a{2z + 1)) 

Aroz{2a + r'^{z^ + z + A)) 

Sr^jz - Ifg^ - 2ro(z - l)(rg(z2 - 2z + 4) + a)gi + 4z{rl{z^ + z + 4) + 3a) 

4roz{2a + r^z"^ + z + 4)) 

4{ z - Ifg^ + Azrljz - l)gi - 2^(rg(^3 + 2^^ + 3^ + 4) + 2a(;^ + 1)) 

2roz{2a + r^iz^ + z + A)) 



C\ — Co 



9\ = 9o- 



Again, in the above expressions one must simply drop the terms highhghted in red 
(or 0" = 0) to get the boundary conditions for the "black brane" case, removing the 
term arising from the S"^. One can easily see that this corresponds to a large tq limit, 
a feature shared with the black hole/black brane solutions in AdS4,. 

2.2 The perturbed solution near r=infinity 

We now turn to the question of the deformation space around the solution given in 
([6]) and (jl]). For this, we will need to consider the term in the potential ae^^^'~^^ to be 
ffist order in e (the perturbative parameter) already. This is because the "background 
solution" is only a solution to the action (1101) with cr = 0. In fact, one can see this 
from the equations of motion written as fl20l) - fl22l) . where the terms multiplied by a 
have more powers of r in the denominator, and so may be neglected in the large r 
limit. Therefore, we take the expansion of the functions 

A{r) = In(r^L) + eAi(r), B{r) = In(rL) + eBi{r) 

C(r) = In + eCi(r), G(r) = In (^^^S^^ ) + eG^{r) (29) 



^We may compare this to the result quoted in [24] and [25] by taking go = ^oCoxT^y 2z{z+i) 

3 

r ^ 

and resolving for cq. This gives cq — , ° In agreement with their results after 

identifying our "cq" is their "go-" 
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and regard the term ae^^^^'^ to be a small correction to the action to find the 
equations of motion 



+ -d [-2{z + l)Ai -{z- l)Bi + -{z - l)Gi + C 

+ 1 ^-l{z^ + z-2)A + l{z^ + z-2)G, + l{z^ + 7z-2)C,^ 

1 

2^ 



= (30) 



-d\B,) + -d -AB, - -{z - l)Gi + Ci 
r \ 2 

+ ^ + z - 2)Ai - h,z^ + z - 2)Gi + h^z^ + z - 6)Ci 

= (31) 

-S^Gi + ^9 - 25i - + 3)Gi + zCi) + = 0. (32) 

In addition, we have the Hamiltonian constraint 

-9 (-2^1 - 2(2 + l)5i - (^ - l)Gi) 
r 

+ ^ ((^' - 3^ + 2)^1 - {z^ - 3z + 2)Gi + {z^ + 3z + 2)Ci) + = 0. (33) 

Note that in the above expressions we have not yet taken the Bi{r) = gauge. We 
turn to this choice shortly. 

We may read the above equations in the following way: The terms with fields in 
them are in fact the perturbed equations associated with action (fTUI) with a = 0, and 
the source terms come from plugging in the 0th order solution into the correction 
term in the action ae^'^^'-'\ These terms are easily identifiable as those being non 
homogeneous in powers of r. 

Using straightforward perturbation theory, we may find the solutions in the Bi = 
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gauge 





+ Ci{z^ + 3z + 2 + {z + 1)7) 



(34) 



+C2 {z^ + 3z + 2- {z + l)-f)r 



2 — zBq + a 



1 



2r^{z^ -2z + 2) 



Ci{r) 





-Co{z - l)r-^-2 ^Ci{z'^ -7z + 6 + {z- 1)7) r 



(35) 
(36) 




1 



2r2(22 -22 + 2) 



(37) 



where we have defined the useful constant 



7 = V9z^ - 20z + 20. 



(38) 



We postpone plotting the exponents of the various powers until section [3] where we 
will discuss how the finite energy condition constrains the constants Co, Ci, C2, Sq. 
The graphs appear in figure [3] of that section. 

The rescaling of time is manifested by shifting Bo, but no such constant seems to 
exist for rescaling Xi (for the a = case). These are in fact gauge equivalent given 
the symmetry of the background r Xr,t — > X~^t,xi^2 ^~^Xi^2, see appendix [Bl 

Above, one must read the powers r~i^^~2 and r~^~5+2 carefully. This is be- 
cause for certain values of z these powers become the same as other modes already 
appearing in the expansion. For example, in the limit that z = 2 the mode r~^~i~^2 
becomes a constant. This constant interferes with the Bo constant term. To get the 
required number of independent integration constants, one must take Bo,Ci to have 
z dependent pieces: the leading orders as 2; — > 2 are engineered to cancel leaving 
the next term in the expansion of lim^^2 ^^ 2"^" 2 (this is a ln(r) piece). Doing so 
consistently will leave the term r~^~S"''2 in Ci(r) to become a constant, even though 
it naively goes to zero when z = 2. Similar considerations appear for the power 
r~2~^~2. In fact, for the z = 2 case, we get modes of the form Bo, ln(r), and 
ln(r)r"*^ using the above considerations [24]. 
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2.3 Finite energy condition. 

Here we explore the finite energy conditions for our backgrounds. We follow the 
discussion in [26] which in turn follows the conventions of [27]. The on-shell Hamil- 
tonian reduces to a boundary integral defined by the intersection of surfaces of con- 
stant "time" and "spatial infinity." For this, one needs to have a (monotonic) time 
coordinate, and for us this is simply t. Further, we need to define a "spatial infinity" 
at a given time slice t; here again, we simply use "r = large constant" to define 
slicing the spacetime near spatial infinity. Further, in the black brane case, one must 
also restrict the integration over xi and X2, and so we must define an energy density 
per unit two volume, rather than an absolute energy. Therefore, we will also need 
the constant x* slices as well, and their associated normal vectors. 
Recall, the metric in the black brane case is 

ds' = -e'^i^Ut' + e^^M (^{dx'Y + {dx'Y) + e^^Wrfr^ (39) 
and so we define the normal vectors to constant t, r, x^ slices as 

nffdf, = e-^'-'^dr (40) 

Next, we define t^d^t = l,f'^d^r = 1, [x^Y d^x^ = 6^^. This then allows us to define 
a lapse function and shift vector for each of these slicings as = [Ntn^ + Nj^) and 
similarly for the other vectors. For us, the shift vectors are trivial [Nj^ = Nl^ = = 
0), and the lapse functions Nt = e'^'^'^\Nj. = e^^'^\ = e^^^'\ These are needed to 
define the measure of the various integrals. We will consider doing a finite integral 
and taking the limits of r_ — > (or to the horizon) and r+ — > oo. In figure [2] we 
give a pictorial representation of the relevant integrals. In [26], it was found that the 
integrated hamiltonian density on a time slice was given bjlfl 

-2 [ {Nt^K - Ntp,,,rY) (41) 

where we have used a generic fj to denote any of the spatial unit vectors (and r* to 
denote which boundary we are talking about: either r = constant, or = constant), 
and Pfj^i, is the momentum conjugate to the time derivative of the metric on the time 

■^In our normalization, we have a coefficient "1" in front of tlie R term in tlie action. Ttierefore, 
to arrive at our normalization, we simply replace the terms IGtt — > l,87r ^ ^ in their expressions. 



12 



Figure 2: The above gives a picture of the cell of spacetime under consideration. r_ 
is the interior surface, r+ is the outer surface with r+ > r_. The red (dashed) line 
gives the surface integral used to define the energy given at the bottom time slice. r_ 
will be taken to approach the horizon (or 0), and r+ will be taken to go to infinity, 
and further Ax* are both held fixed, with their product Ax^ x Ax^ = V2. 



slice. For us, this second term is not present because all shift vectors N'^^^.^ 2} = 0- 
The term is the extrinsic curvature of the spatial boundary slice in the constant 
time slice t = constant. We therefore have 6 integrals to perform, for the 2 limiting 
values of (r, x^,x^). This gives 



I r=r+ 
I r=r— 



-2 J dx2dr [Nt {g^" + nf^ - n^,K) V,n,,] ifCllt (42) 
-2 J dx.dr [Nt {g^- + nf< - n^<) V .n^,] 



^2- 



where the bounds of any x* integration are (x*~,x*"'" and that of any r integration 
are (r~, r"*"). The last two terms above are zero because the integrand is independent 
of x^ and x^. Further, the first integrand is independent of x^ and x^ and so the 
integral is just V2. Recall that x* are unitless, so at the end we will need to take 
some factors of L along with this unitless V2 to honestly get an energy density. This 
then simplifies to 

= -2[e^Me-^W(9,(e2^M))]|;z;+. (43) 

In the above equation the combination e~'-'dr is r-diffeomorphism invariant, and so 
the above is well defined for any redefinition of the r coordinate. In the pure Lifshitz 
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background the limit r_ — > gives a zero answer. Further, we note that for the black 
hole backgrounds we are dealing with, 0, e~'^^^^ — > 0, and e^^^''^ a constant 

when r ^ vq. This means that in both cases, when one takes r„ to its limiting value, 
this boundary will not contribute. Therefore, only the boundary at r = r+ — > oo will 
contribute. To render this finite, we will use a background subtraction technique. 
Keeping in mind that the zeroth order Lifshitz solution is what we compare to, and 
working in the -Bi = gauge, we will find that the r~^^^ mode is the mode that 
has finite energy. Any function that does not fall off at least this fast will give 
an infinite energy density contribution, even if the background asymptotes to the 
Lifshitz solution. We will discuss this further in the next section. 

One may repeat the above calculation for a = 1. One need not consider taking 
finite two volume in this case because the S"^ is compact. Therefore, the analogs of 
the last four (out of six) terms of ( H2l) are not present. Further, in equation ( l43l) one 
replaces S^^ — g^z- However, it is still true that F^^- = —^d (e^*^^^) g{s^)ij- Given this, 
the same formula above applies, with the simple replacement V2 —>■ 47r to account 
for the finite volume of the unit S*^ . 

Again, we use the B{r) = In(Lr), Bi{r) = gauge. In such a situation, the above 
energy reads 

p = lim -4 re^(")e-^("Vl (44) 

for all backgrounds. 



3 Numeric integration, results, and discussion. 

First, we will consider the finite energy condition applied to the large r region of 
our backgrounds. We subtract a reference background to obtain a relative energy 
density 

p,ei = lim -4 \r (e^^^^e-^^^^ - e'^-^^'^e-^-'^'^)] . (45) 

We require both backgrounds to be asymptotic to the Lifshitz fixed point. Therefore, 
we expect that for r ^ cxd the perturbation theory developed in section 12.21 applies. 
This gives 



Prei = lim -4 \r (e^^'-^e-^^'') - e^-^^^^e-^-f^'"))] 

r_|_— >cxD 

lim -4 _ _ (A,_,,f(r) - Ci,,ef(r)))] . (46) 
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a (z) 



1.0 1.5 2.0 2.5 3.0 

z 



-5 

-6 

Figure 3: Above we graph the exponents of the modes at infinity r"'*^^^ as a function 
of z. The (sohd) black fine {—z — 2) represents the proper "energy" mode. The (long- 
dashed) red line (—2) represents the inhomogeneous mode induced by the S"^. The 
line (—2), and all its non-linear descendants, will be universal for backgrounds with 
the sphere. All other curves that are above the black line {—z — 2) represent infinite 
energy modes, and curves below this line represent modes with contribution to the 
energy of the background. The top (short-dashed) magenta curve is — 1 + 
(recall this mode has coefficient Ci), which always represents an infinite energy mode. 
The bottom (dash-dot) blue curve is — 1 — j (recall this mode has coefficient 

C2). This curve represents a finite energy mode for z > 2 and an infinite energy mode 
for 1 < 2; < 2: the value z = 2 gives the curve where this curve crosses the line —z — 2. 



We plot the exponents r"*-^) in the modes of flM|) - fl37|) in figure [31 We first discuss 
the (7 = case, where the term is not present, and in fact we may use the regular 
Lifshitz background as our reference (yliref(r) = Ci^i.cf('") = 0). In such a case we 
require that the coefficient of r~t~^+2 be zero, Ci = 0, so that the limit (H6l) does not 
diverge. From figure [3l we see that after setting Ci = for z > 2, the next to leading 
mode at infinity is the r^^^^ mode. This precisely cancels the r^"*"^ term in (1461) . 
so that the r — > 00 limit is finite. Hence, for z > 2 we need to cancel one term at 
infinity: we must set Ci = 0. Now recall that we have two parameters at the horizon, 
and therefore we expect a one parameter family of solutions with finite energy. This 
single parameter is most naturally taken to be the position of the horizon. 
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However, in the case that (1 < z < 2), there is an additional mode we must 
remove, namely r^5~^^2. This gives two conditions at infinity. Recall that we 
always have two constants at the location of the horizon. Hence, on general grounds 
we expect that there is at most a discrete set of solutions that have a regular horizon 
and have finite energy density in the regime [1 < z < 2). However, this is actually 
not the case for a = 0. 

For the a = black branes, the scaling symmetry t X^t, x' — > Xx\ r X^^r is 
exact for the background: only the presence of the horizon breaks this. This implies 
that in the a = case, if there exists one finite energy black brane (i.e. one finds 
one combination of go and rg that gives a finite energy density black brane), with 
horizon position Tq, then a continuous one parameter family of solutions must exist. 
This is because we may apply the scaling symmetry and find a new finite energy 
black brane, continuously scaling the location of the horizon Tq — * Xtq. One can see 
this most easily by considering the correction functions in the metric and form field. 
Assume we have a finite energy density solution 

ds^ = (^-r^^e'^'^^'''-'^'''^df + r'idxj + dxj) + e^^^^^"'^"'^^)^^ 
z \ 

After performing the rescaling, the solution reads 

Now, if e^^^^^°''''°''^^ has a linear zero at r = rn, then e^'^^'^s'o.fo.A ^r) j^^^g ^ linear zero 



at r = Xtq (and so forth for the other functions). Further, no new terms in the 
expansion at infinity have been introduced. Hence, the scaled solution is also a finite 
energy density black brane with respect to the uncorrected Lifshitz background. 

Thus, this scaling symmetry allows us to conclude the following result for the 
black brane case: either there is a 1 parameter class of black branes (i.e. Ci = 
implies C2 = for physical input on the horizon 5'o;''"o); or there are none at all. It 
would be interesting to see if there are any by shooting in from infinity and seeing 



"^we include the ~ in these bounds because we get logarithmic terms at the end points: these 
still have logarithmically divergent energy density 
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if a regular horizon develops, given a finite energy perturbation at infinity. If one 
does not see a regular horizon, then one may conclude that there are no finite energy 
density black brane solutions, and conversely if there is one, then there is a continuous 
family. We will see later that for the z > 2 case a continuous class of black branes 
exist. However, applying the above scaling argument to the 1 < z < 2 case, we see 
that the naive expectation that there are a discrete set of black branes is incorrect; 
actually there are either none, or a continuous class of them. 

At exactly z = 1, the counting of constants at infinity is different because one- 
form field is set to zero. In such a situation one may count 2 dynamical functions. As 
we saw in the perturbative expansion, the Hamiltonian constraint and gauge fixing 
remove two additional constants, leaving two constants at infinity. One of these 
must be associated with time diffeomorphisms, and so we arrive at only having one 
meaningful constant at infinity. Setting z = 1, we also only have one constant at 
the location of the horizon. Therefore, if the mode found at infinity is finite energy, 
then one expects a one parameter family of solutions. If one wishes to turn on the 
one-form field, one would instead count 2 independent modes at the horizon and 3 
at infinity. This counting is of course valid for both cr = 1 or cr = 0. 

Much of the above discussion is directly applicable to the a = 1 case. For z > 2 
we again see that we must cancel the top mode Ci = 0. We also have the term, 
and all of its higher order corrections to deal with. However, the comes with a 
known coefficient, and so do all of its "descendent" modes. So, we get r~^,r~^, etc. 
terms, all with known coefficients, which are the same for all solutions. Therefore, 
we expect these modes cancel against any given reference a = 1 background that 
asymptotes to the Lifshitz solution (see for example the exact solution in appendix 
|A]or the ^^5*4 black hole/black brane solutions). Again, for z > 2, this cancels all 
infinite energy modes. Hence, for 2; > 2 we expect to obtain a one parameter family 
of finite energy solutions. 

When 1 < 2; < 2 we again find that we must set both Ci = and C2 = 0, 
constraining the two free parameters at the horizon. Because of this, we expect to 
generically find at most a discrete set of solutions. This would have been true even 
if we took arbitrary Ci,C2 as the background solution: fixing Ci = 0,C2 = is just 
a particular case of this. Let us consider looking for backgrounds with fixed values 
Ci^ref and C2.ref- There are still two conditions at infinity and so we still need to 
tune both go and tq so that both Ci and C2 match their reference values. Therefore 
we would still expect at most a discrete set of go, tq that satisfy these conditions. 
Other differences between the z > 2 and 1 < z < 2 were already found in [28]. The 
presence of the sphere in this case breaks the scaling symmetry. The lack of a scaling 
argument does not allow us to make as strong a statement as for the cr = case. So, 
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generically we expect a discrete set of solutions forcT = l, 1 < z < 2, analogous to 
the discrete set of Lifshitz stars found in [24] . 

Given the above considerations, the program for finding black holes/branes is as 
follows. We have two free parameters at the horizon, qq and tq. We consider fixing 
tq and scanning through possible values of go until we find a solution that at infinity 
admits the expansion ( l34l) -( !37l) . We further scan through go until the coefficient 
Ci = for our solutions, as we expect this mode to either give divergent behavior 
(for 2; > 2) or have infinite energy (for any value of z). This is sufficient to find finite 
energy (density) black holes/branes for the z > 2 cases. For the 1 < 2; < 2 cases, 
we would also have to scan through rg until C2 was zero as well. However, we will 
be unable to do so, as we can only set Ci numerically close to zero, and not actually 
zero. We will not be sensitive to the power in front of C2 because the non linear 
corrections to the linearized equations. 

One could in principle do a higher order in perturbation theory calculation, and 
go to sufficiently high order that one could be sensitive to C2. We may get a ballpark 
estimate of how high in perturbation theory one would need to go before one would be 
numerically sensitive to these modes. Let us take 2; = 3/2 as an example. In this case 
— | — 1 + 1 ^ —0.1492 and —| — 1 — I ~ —3.3508. Generically, this term is present and 
we can only numerically set its coefficient Ci to some small number. Hence, to find 
out what the value of C2 is, we would need to remove r~2~i+2 and all its subsequent 
descendent modes arising from nonlinearity in the system. These descendent modes 
will typically be integral powers of r~S~^+2, and so we ask the question, for what 
integer n is n (— | — 1 + 2) (~f ~ -'- ~ 2) that we are sensitive to the power 
multiplying C2. For the case of 2; = 3/2, the answer is n = 23, and so one would have 
to go to at least 22""' order in perturbation theory separate out C2! For this reason, 
we do not attempt to do this here. In fact closer to z = 1 and z = 2 this problem 
becomes worse because ~f ~ 1 + 2 approaches 0. 

In what follows, we will simply set Ci to zero, knowing this is sufficient for the 
z > 2 cases and necessary for the 1 < z < 2 cases. We will also restrict to the 
(7 = 1 case: the a = solutions can be obtained from the plots below by simply 
looking at large tq: in this situation the term e^^'^^^'~'^^^ potential term always remains 
perturbative, and so is arbitrarily close to setting a = 0. Setting Ci = will furnish 
go in terms of tq and so we can find the temperature via 

T = (49) 

4vrco 

where oq is chosen so that e"^^^^-* asymptotes to 1 rather than an arbitrary constant 
and Co is given in terms of go and tq in ( l25l) . Of course this is a "unitless" temperature: 
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all units are restored with L. 

Again, we stress that we do not have the numeric precision to be able to set the 
coefficient C2 = for the 1 < z < 2 cases, as explained above. Because of this extra 
constraint, we expect to get a discrete set of black holes. Therefore, for 1 < z < 2 
the temperature T(ro) plots are to be regarded as the curve along which a discrete 
set of points represent good finite energy black holes; for 2; > 2 the curves represent 
a one parameter family of black holes. 
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Figure 4: Graphs of log^g (^(^0)) for ^ = 4, 3, 2.01 colored red (solid), green (dash- 
dot) and blue (dashed) respectively for the a = \ case. Note the feature developing 
for smaller values of z. 

There is a possible exception to this: if for some reason setting Ci = also gives 
C2 = given the physical boundary conditions at the horizon, then these curves are 
actually the temperatures of a continuous set of black holes. Although we feel that 
such a situation would be non generic, we leave it open as a possibility. In figure H] 
we plot the T(ro) for z > 2 cases z = 4,3, 2.01. 

In figure m a flattening of the functions occurs at logiQ(ro) ~ —0.5 for smaller 
z. One may wonder whether these curves develop a negative slope, and so signal 
some sort of thermodynamic instability. In figure [5a] we plot the T(ro) curve for 
2; = 1.1 and compare this to the pure AdS black hole, and indeed we find that near 
log^o('^o) ~ —0.25 the slope becomes negative. If this represents a continuous set of 
black holes this would signal a thermodynamic instability: if however it is just the 
curve along which a discrete set of black holes may exist, then the interpretation is 
not obvious. 
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1.1, z = 1 




(b) z = 1.74,1.761,1.78 



Figure 5: On the left is graphs of logj^Q (T(ro)) for 2; = 1.1 colored red (dash-dot), and 
the temperature for z = 1 the normal AdS black hole with temperature T = 
in green (sohd). On the right, we plot log^Q (T(ro)) for z = 1.74, 1.761, 1.78 in blue 
(dashed), green (solid), and red (dashdot) respectively. All graphs are for the a = 1 
case. 



We also find the critical value of z where a negative slope of T(ro) is first devel- 
oped; we find this to be 2; ^ 1.761. This is shown in figure [5bl One further note 
is in order. Although T(ro) develops a negative slope for certain values of tq, for 
sufficiently small tq it seems to become positive again. If we take the negative slope 
as a sign of a thermodynamic instability, it does not have the runaway behavior of 
the pure AdS black holes. It is eventually "caught" from shrinking further after it 
becomes sufficiently small. This is, of course, assuming that we have a continuous 
set of black holes. 

For completeness, we also show several plots of the functions e^^^^^^ e*"^^*"^ and 
gGi(r) figure O One should notice the universal behavior for large tq; the functions 
remain monotonic, and that for small tq, that they develop an extra feature. This is 
similar to the behaviors of small and large black holes in AdS, which are also graphed 
in figure [6l Again, the large tq behavior of all the graphs is expected to be identical 
with the cr = CcLSG, clS the term in the potential e^^*" always remains perturbative. 
So we do not consider this a separate case: simply a limit of the solutions given here. 

We conclude with some open problems. The most striking is the question of 
whether for 1 < z < 2 finite energy black holes form a continuous family, or a discrete 
set. This seems to be a difficult question to answer, even numerically, given the order 
in perturbation theory one would need to attain (discussed above). Perhaps the most 
efficient way to address this problem is shooting in from finite energy conditions at 
infinity, and determining whether one develops a regular horizon. 
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(c) To = 0.1, pure AdS 



(d) tq = 5, pure AdS 



Figure 6: Graphs of the black hole correction functions e^^^^\ e'"^^^\ e'-^^^'^\ colored 
red (short-dashed), blue (dash-dot) and green (long-dashed) respectively, for the 
cr = 1 case. The (solid) black flat line is to guide the eye to the asymptotic value 
of 1 for all the functions. Note that the large black holes remain monotonic, while 
the small black holes have a feature at several times tq. Further, it takes many more 
times ro for the solution to asymptote to 1 in the "small" case. This behavior is 
expected from the black hole in AdS, also graphed above for comparison. The black 
brane a = case is given by the large ro case above. 



One could also consider other "finite energy" conditions that one might be able to 
impose on the above solutions. One possible candidate is the program of using local 
counter terms at the boundary [29,30], to remove the infinite energy contributions for 
black hole backgrounds [31]. Such a program might give some connection between 
the modes multiplied by Ci and €2- 

Of further interest would be to embed the above action and/or solutions into 
a supergravity. This could give some clue as to how these solutions are embedded 
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into string theory, but also may open up new branches of solutions by introducing 
new matter content, and further one could look for possible supersymmetric versions 
of these backgrounds. Particularly interesting is the construction of [32] where non 
relativistic spacetimes were realized as cosets of the Schrodinger algebra. This raises 
the interesting question of whether these solutions can be constructed as the bosonic 
part of some supercoset models. We look forward to addressing these issues in the 
future. 

note added: 

After releasing our work, the program of constructing local boundary counter terms 
has been performed [33]. There, the authors show that when the largest mode at 
infinity is set to zero (Ci = 0) the energy is finite for any z. This leads us to conclude 
that the above black holes for 1 < 2; < 2 are in fact a continuous set. 
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A An exact solution. 

In this section we will present an exact solution (we really mean an: there are no free 
parameters). In fact, we are able to find 2 solutions, one with a naked singularity. 
Since both of these solutions follow from similar guesswork, we will display them 
simultaneously. If we start with equation ( l23l) and require that e^^'<^^ e^^'<^^ = 1 (or 
actually, any constant, but this may be removed with rescaling time), we in fact get 
the algebraic relation e*-^^ = ±e^^^. This we plug into ( 122|) and get a ear equation 
in e^^^. It reads 

- r^dde^^'^'^ - r{z + 3)de'^'^'-^ - 2ze^'^'^'^ + 2z (50) 
and has the solution 

e'Mr) = M,r-^ + ^ + 1 (51) 
which converges to 1 for large r. We then plug this into either (j20|) or (12T]) and find 

-2z + {4z - 4z^ + 2z^)M2 + ^ (52) 
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This can obviously not be solved for generic z. However, for z = 2 and for 2; = 4 we 
do find solutions. 




e-2C.(.) ^ ^GAr) = ^2Mr) = ^ ^ " ^ (53) 

^ . . 1 3 if ^ = 4 ^ > 



400r4 



Clearly the 2 = 4 case has a horizon. The z = 2 case, however, does not (a recent 
paper [25] also has this solution). One may be concerned about the fact that the 
two modes l/r^ and are the same in this case, and one must really resolve the 
equations and get a logarithmic piece. However, one can show that the coefficient to 
the logarithmic piece must be zero, and so one simply arrives at the l/(2r2) + l type 
of solution. One can show that the z = 2 solution has a curvature singularity at r = 
and that this at finite (spacelike) geodesic distance from any point in the space. We 
therefore consider this an unphysical solution. However, the z = 4 solution is a good 
solution, with finite Ricci tensor, and finite Riemann square for all r > ro = 
where this value of tq is the position of the horizon. Both the Ricci tensor Riemann 
square tensors only become infinite at r = 0. 

For completeness, we write the metric and one-form explicitly here. 



ds^ = -r-s ( 1 + -_ - — _ 1 de + r'd^i' + 



dr"^ 



lOr^ 400rV + ^ ' ^) 

^ =Lr\ll{l + ^,~-^^dt (54) 



2 V 10r2 400r4^ 
where dVL^ is the line element on the unit two sphere 



B Gauge invariance 

In the previous sections, we have gauge fixed by taking Bi{r) = 0. Here, we write 
down the linearized gauge transformations that will allow us to switch to other gauges 
in perturbation theory (used near r = oo). The transformation 

Ai(r) Ai(r) + -5(r) 

r 

5i(r) ^Bi{r) + -6{r) (55) 
r 

Ci(r) ^Ciir)--6{r) + drS{r) 
r 

Gi(r) -.Gi(r) + -5(r) 

r 
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corresponds to infinitesimal coordinate transformations r ^ r + e(5(r). One can 
see tliat sucli a sliift leaves the first order equations (near r = oo) unchanged by 
explicitly plugging in the above shift. 
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